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This talk summarizes recent QCD results from HERA, the Tevatron Collider and Tevatron fixed target experiments. 
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1 Introduction 

As implied by the title of this talk, the topics to be dis- 
cussed cover a very wide range, encompassing QCD re- 
sults from the Tevatron pp collider, the HERA ep collider 
and the Tevatron fixed target experiments. In my talk, 

1 will not try for a totally comprehensive review, but 
instead will discuss some of the important experimen- 
tal and phenomenological developments in perturbative 
QCD since the last Rochester conference. I will not be 
covering deep inelastic scattering results per se, which 
will be discussed in the summary talk of Tony Doyle; 
instead I will present results which emphasize the de- 
tails of the hadronic final state in DIS and photopro- 
duction events. Similarly, I will not discuss diffraction 
which will be the subject of the review talk by Martin 
Erdmann. Many of the recent developments in the the- 
ory/phenomenology world, along with a discusion of our 
current understanding of a s , will be contained in the talk 
of Yuri Dokshitzer. 

The main theme of my talk will be the success with 
which perturbative QCD has been applied to the data 
from Fermilab and from HERA. There are enough mys- 
teries left, however, to make life interesting (and to pro- 
voke the need for larger data samples), with several of 
the mysteries involving the remaining uncertainties in 
the gluon distribution. DGLAP-based perturbative QCD 
predictions remain very successful and the search for con- 
vincing evidence for BFKL effects continues. 

2 Tevatron Collider 

In the Tevatron collider, 900 GeV protons collide with 
900 GeV antiprotons leading to a center-of-mass energy 
of 1.8 TeV, the highest energy currently accessible. The 
Tevatron collider completed a very successful Run 1 in 
1996 with each experiment (CDF and DO) accumulating 
on the order of 100 pb' 1 of data. Most analyses have 
been published or are nearing publication. 

2.1 Inclusive Jet Production at the Tevatron 

The inclusive jet cross section in the central rapidity re- 
gion has been measured by both the CDF and DO exper- 



iments at a center of mass energy of 1.8 TeV. Jets are 
defined using an iterative fixecLcone algorithm with a ra- 
dius R (yj An 2 + A(f> 2 ) of 0.7. El The measurement spans 
the transverse energy range from 15 GeV/c to the order 
of 500 GeV/c; in this range the jet cross section drops by 
over 9 orders of magnitude. The highest Et jet events 
probe the smallest distance scales (10 -17 cm) currently 
accessible. Any new physics that might exist at these 
distance scales, such as compositeness, might manifest 
itself in the jet cross section measurement. 

The jet cross sections from both experiments are cor- 
rected for detector measurement and resolution effects 
and are compared to next-to-leading order (NLO) QCD 
calculations. opij The theoretical uncertainties in inclu- 
sive jet production are greatly reduced at NLO as com- 
pared to leading order. The. two programs that are cur- 
rently in use are JETRAD L3 and EKS. I JETRAD gen- 
erates the NLO inclusive jet cross section by a Monte 
Carlo phase space slicing technique while EKS is an an- 
alytical calculation. Both programs are implementations 
of the same matrix elements^ and give essentially equiv- 
alent results when the same cuts/conditions are applied. 
At NLO, the sensitivity of the jet cross section to the 
renormalization/factorization scale is reduced, but still 
present. The value of this scale should be proportional 
to the hardness of the scatter. It is convenient to set the 
rcnormalization and factorization scale equal to a multi- 
ple of the Et of the measured jet. One can also set these 
scales to a multiple of the maximum of all jets in the 
event (Exmax), at the cost of introducing another vari- 
able into the prescription. Typically, ETjet/2 is used in 
the EKS program while E Tmax /2 is used in JETRAD. | 
Exjet/2 may be a more "natural" choice for an inclusive 
jet cross section, but E Tmax /2 is also acceptable. The 
use of E Tmax /2 rather than E T j et /2 leads to a reduction 
in the jet cross section of 7% at Et — 50 GeV/c decreas- 



Tn the JETRAD program, each Monte Carlo event has 2 or 3 
partons in the final state, leading to a possibility of either 2 or 
3 jets. The Et of each individual jet is not known until a jet 
clustering algorithm has been applied (the two lower Et partons 
may be clustered if they are close together). The only scale known 
unambiguously at the time of the event generation is the Et of 
the most energetic parton (ETmax)- A version of JETRAD also 
exists in which a second pass is made through the generated events, 
allowing the use of the scale ETjet ■ 
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Figure 1: The inclusive jet cross section from CDF for Run 1A and 
Run IB compared to the NLO QCD prediction using the CTEQ3M 
parton distributions. 



Figure 2: The CDF inclusive jet cross section from Run IB com- 
pared to NLO QCD predictions using the CTEQ4M and CTEQ4HJ 
parton distributions. 



ing to < 1% at E T = 100 GeV/c. The effect on the NLO 
inclusive jet cross section of variations in the renormaliza- 
tion/factorization scale, the value of R se p (the minimum 
separation of the two partons for them to be considered 
as two separate jets), and the choice of parton distribu- 
tion functions (pdf's) is investigated in more detail in 
Ref. 6. 

The picture that has been adopted for jet production 
is that the final state consists of 2 or 3 partons from the 
hard scatter accompanied by an underlying event due to 
the collision of the proton and antiproton remnants. The 
underlying event is taken to be identical to that observed 
in minimum bias events and its contribution to the jet 
cone energy is subtracted before any comparisons to the- 
ory. This picture has been successful but may be incom- 
plete; there may be additional contributions to the un- 
derlying event from double parton scattering and higher 
order radiation effects that may not be included in the 
subtracted minimum bias level and may not he correctly 
described by the NLO QCD calculations. oBEl A uniform 
extra contribution to the jet energy might manifest it- 
self as a jet profile broader in experiment than in theory. 
Such an effect has already been observed by both CDF 
and DO. The main impact of any underestimate of the 
underlying energy level in jet events would be on lower 
energy jets. 

The CDF collaboration has previously published the 



inclusive jet cross section from Run 1A (19.1 pb^ 1 ) for 
0.1 < |r/| < 0.7. E£l A linear comparison of (Data- 
Theory) /Theory is shown in Figure [l] along with the pre- 
liminary results from Run 1B(87 pb^ 1 ). Good agreement 
with the NLO prediction is observed except at the high- 
est values of transverse energy. The excess is inconsistent 
with the (highly correlated) systematic error of the mea- 
surement and cannot be explained by a different choice of 
renormalization and/or factorization scale, or by a differ- 
ent choice of conventional parton distribution function. 

The CTEQ collaboration has performed a global pdf 
fit using the Run 1A jet data from CDF, giving a large 
emphasis in the fit to the high Et data. The resulting fit 
(CTEQ4HJ) contains a gluon distribution substantially 
greater (by aJactor of 2 at x=0.5) than that in conven- 
tional pdf's. EiJ The larger gluon distribution at high x 
leads to a greater cross section at high Et (20% at 450 
GeV/c). The increased jet cross section prediction does 
not pass directly through the center of the high Et data 
points, but does pass through the bottom of the error 
bars. The CDF jet cross section from Run IB is shown 
in Figure|| compared to the NLO QCD predictions using 
the CTEQ4M and CTEQ4HJ pdf's. 

For high Et jet production, the dominant subprocess 
is qq scattering. The gq subprocess comprises only 20% 
of the cross section and the gg subprocess contribution 
is minimal. For this reason, an increase in the gluon 
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Figure 3: The ratio of gluon distributions consistent with the DIS 
and DY data sets to the gluon distribution from CTEQ4M. The 
gluon distribution from CTEQ4HJ is also shown for comparison. 



distribution of a factor of 2 leads to only a 20% increase in 
the jet cross section. The quark distributions are also free 
to change in this fit, but are tightly constrained by the 
highly precise deep inelastic scattering (DIS) and Drell- 
Yan (DY) data at these x values. 

It may seem suprising that the gluon distribution has 
this degree of flexibility. A recent CTEQ paper explored 
the uncertainty in the gluon distribution by performing 
a gluon parameter scan,_utilizmg the DIS and DY data 
used in the CTEQ4 fit.B The resulting pdf's were ex- 
cluded if there were any clear conflicts with any of the 
data sets. The pdf's that remain are shown in Figure || 
and indicate that the gluon is tightly constained at lower 
x. DIS and DY data provide little constraint, however, on 
the high x gluon distribution and this is demonstrated in 
the much wider variation observed in the large x region. 
(The CTEQ4HJ gluon distribution is also indicated for 
comparison purposes.) In previous pdf's this constraint 
has been provided either by fixed target direct photon 
data □ and/or jet data from the Tevatron collider. 
Due to evolution effects, the gluon distribution is more 
tightly constrained (for x < 0.2) at high Q 2 than at low 
Q 2 - 

DO has presented at this conference a measurement 
of the inclusive jet cross section for < 0.5 (and 
0.1 < \rj\ < 0.7 for direct comparison .to-CDF) with 
a substantially reduced systematic error. Iia€3 The cross 
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Figure 4: A comparison of the DO inclusive jet cross section and 
NLO QCD (JETRAD) predictions obtained using three different 
parton distribution sets. The band represents the total experimen- 
tal uncertainty. 
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Figure 5: A comparison of the DO inclusive jet cross section and 
a NLO QCD (JETRAD) prediction using the CTEQ4HJ parton 
distribution set. 

section for \rj\ < 0.5 is compared to NLO QCD predic- 
tions with several pdf's in Figure [|. Good agreement is 
observed with perhaps some sign of an excess at moder- 
ate to high E T when the CTEQ3M and CTEQ4M pdf's 
are used. The data is uniformly larger than the predic- 
tion using the MRST pdf. This latter deviation is due 
to the substantially weaker MRST gluon distribution at 
high x (see Section 3.1). 

Taking into account the correlations in the system- 
atic errors, reasonable \ 2 values are obtained for all 3 
pdf's. The MRST prediction agrees in shape with the DO 
data, although not in normalization. The best \ 2 agree- 
ment (nominally better than MRST) is obtained with the 
CTEQ4HJ pdf (shown in Figure | for O.K |r?| < 0.7). 

A reanalysis of data from SLAC and NMC, taking 
into account nuclear binding effects in the deuteron, pre- 
dicts a larger d anark distribution at high x than found 
in modern pdf's. tZI One of the consequences of this larger 
high x d quark distribution would be an enhanced high 
Et jet cross section (by about 10% at the highest Et 
values) . 
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Figure 6: The top plot shows the normalized comparison of the DO 
inclusive jet cross section to NLO QCD predictions (EKS) along 
with a similar comparison of the Run 1A CDF inclusive jet cross 
section. The CTEQ3M parton distribution set is used. The mid- 
dle plot indicates the difference between the DO inclusive jet cross 
section and a smooth fit through the Run 1A CDF inclusive jet 
cross section, normalized to the latter. The band represents the 
uncertainty on the DO data. A comparison of the CDF and DO jet 
systematic errors is shown at the bottom. 



A direct comparison of the CDF Run 1A and DO jet 
cross sections is shown in Figure ^. Aside from a nor- 
malization shift of 5%, the two experiments obtain very 
similar jet cross sections, with the most noticeable differ- 
ence being in the last two data points. A normalization 
shift of this magnitude is to be expected since CDF and 
DO use different values for the total inelastic pp cross 
section. DO has calculated a low probability for their jet 
data to agree with the physics curve drawn through the 
CDF Run IA data, but if the correlated errors for both 
experiments are taken into account, there is agreement 
between the two at the 46% level of probability. 

2.2 Dijet Cross Sections at the Tevatron 

Both experiments have reported results for the dijet mass 
cross section and observe agreement with each other and 
with NLO predictions, albeit with an excess at high dijet 
mass consistent with that predicted by the CTEQ4HJ 
pdf , as shown in Figure Il3t3 

DO has used the dijet mass cross section ratio of \rj\ < 
0.5 over 0.5 < \rj\ < 1.0 to test for compositeness. At the 
95% CL lipit, a compositeness scale of 2.4 TeV can be 
excluded. E3 Previous measurements of the dijet angular 
distribution provided an exclusion of A < 2.0 TeV. 

CDF has presented a measurement of the differential 
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Figure 7: A comparison of the CDF(preliminary) and DO 19 di- 
jet mass cross sections to the NLO QCD predictions using the 
CTEQ4M, CTEQ4HJ and MRST parton distribution sets. 



dijet cross section in which one jet (the trigger jet) is re- 
quired to be central (0.1 < \r\\ < 0.7), while the other j£t 
(the probe jet) can have any rapidity value up to 3.0.113 
The differential dijet cross section is then plotted versus 
the transverse energy of the trigger jet, for the 4 different 
probe jet rapidity intervals. The measurement presented 
in this way takes best advantage of the better jet Et reso- 
lution CDF has in the central rapidity region. This mea- 
surement also directly probes higher x values than the 
inclusive jet cross section. The dijet cross section is sen- 
sitive to the high x gluon distribution and anything un- 
usual that may occur at high x and Q 2 . This can be seen 
from Figure || where i{= 2E T cosh 2 (r]*)(l - tanh(i]* j) is 
plotted versus x max for the dijet cross section bins. The 
box in the upper right-hand corner indicates the region 
of phase space where a possible excess at HERA has been 
probed. 

A comparison of the measured dijet cross sections to 
predictions using the CTEQ4M, CTEQ4HJ and MRST 
pdf 's is shown in Figure |^. An excess is observed at high 
Et (corresponding to high x) for each of the probe jet 
rapidity bins. The size of the excess is decreased when the 
CTEQ4HJ pdf is used. A detailed conclusion may wait, 
though, until a more detailed study of the adequacy of 
the NLO predictions is carried out for the high Et, high rj 



b For those of you who still do not get the Prince joke, please send 
email to me at the address given. 
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Figure 8: A plot of the t vs x ma x reach for the CDF differential 
dijet analysis. The box in the upper right hand corner indicates 
the kinematic region where an possible excess at HERA has been 
probed. 



region. (At very high x, multijet final states are common 
and NLO phase space may not be adequate.) The MRST 
predictions are uniformly below the data at all Et values 
due to the weaker gluon distribution discussed previously. 

Whenever two jets of roughly equal Et values are 
separated by a large rapidity interval, the emission of 
gluons in the rapidity region between the two jets gen- 
erates logarithmic contributions [(a s \n(s/p T )) n ] to the 
dijet cross section_.which need to be resummed using the 
BFKL equation. E3 Naive BFKL predictions lead to be- 
havior that differs dramatically from that obtained from 
fixed-order perturbative predictions. The imposition of 
kinematic constraints on the calculations, however, su- 
presses the BFKL-like behavior at the Tevatron and re- 
duces the chances of unambiguous observation of BFKL 
signatures. The kinematic environment at the LHC will 
be more favorable for observation of BFKL-like behav- 
ior, m 



2.3 Jet Production at 630 Ge V and the Xt Scaling Ratio 

CDF and DO have both presented measurements of the 
inclusive i£jt cross section at 630 GeV, and of the Xt scal- 
ing ratio. Ej A comparison of the cross sections of the two 
experiments to NLO predictions is shown in Figure [l~0| . 
For both experiments, deviations from the NLO predic- 
tions are observed for jet Et values below 90-100 GeV/c. 
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Figure 9: A comparison of the CDF differential dijet cross section to 
NLO QCD (JETRAD) predictions using the CTEQ4M, CTEQ4HJ 
and MRST parton distribution sets. The systematic error band is 
indicated at the bottom for each probe jet rapidity interval. 



When the scaled cross section ratio (1 / \2n)E T da / 'dEx 
for 630/1800) is plotted versus xt {—2Et/ Vs), many of 
the systematic uncertainties cancel for both experiment 
and theory. A naive parton model prediction would give 
a value for the ratio of 1.0; QCD effects change the pre- 
diction to closer to 2 (with some dependence on xt as ob- 
served). Both experiments have measured an xt scaling 
ratio lower than the theoretical prediction for xt < 0.3. 
(See Figure O.) A similar discrepancy was observed for 
an earlicr_comparison of 546 GeV jet data to 1800 GeV 
jet data.E] The reason for the discrepancy is still under 
theoretical investigation and may be due to a combina- 
tion of effects (underlying event subtraction, initial state 
kr and additional nqa-perturbative jet fragmentation ef- 
fects ("splashout"))J 

2.4 W + Jet(s) Production at the Tevatron 

Another ratio of observables with reduced systematic er- 
rors is the ratio of W + jet(s) production to W produc- 
tion. This measurement is naively sensitive to the value 
of a s and, in fact, was originally proposed as a means 
of measuring a s . DO has reported an exclusive measure- 
ment of W+l jet production to W + jet production at 
several recent conferences. The jets were measured with 
the standard DO jet algorithm using a cone radius of 0.7. 
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Figure 10: A comparison of the CDF and DO inclusive jet cross 
sections at 630 GeV to NLO QCD predictions using the MRSA' 
parton distribution set. The shaded band indicates the DO system- 
atic error. 



The result was in serious disagreement with the NLO 
QCD predictions and no choice of scale or pdf provided 
any significant improvement. 

CDF has measured the ratio of W+>1 jet to inclxLj 
sive W production (Rw) using a jet .cone radius of 0.4l3 
and (new for this conference) 0.7. Ej A comparison of 
both cone size results to NLO QCD predictions Ej in 
Figure |lj indicates good agreement. In Figure |l3|, the 
CDF data is compared to NLO predictions using a va- 
riety of pdf's corresponding to different a s values. The 
theoretical predictions have a surprisingly small depen- 
dence on the value of a s . The experimental ratio of R±q 
with 0.7 cones to Rio with 0.4 cones is shown in Figure |l4] 
along with the theoretical predictions. The experimental 
ratio is larger than the theoretical one, which is another 
indication that jets at the Tevatron are broader than the 
theoretical predictions. A subtraction of an extra under- 
lying event energy contribution would improve the agree- 
ment of the experimental jet shape with theory, and thus 
improve the agreement of the experimental i?io(0.7/0.4) 
ratio with the theoretical one. 



3 Tevatron Fixed Target 



Figure 11: A comparison of the scaled inclusive jet cross sections 
(630/1800) for CDF and DO to NLO QCD predictions. The shaded 
band gives the DO systematic errors. 

3.1 Direct Photon Production 

Direct photon production has long been viewed as an 
ideal vehicle for measuring the gluon distribution in the 
proton. Ej The quark-gluon Compton scattering subpro- 
cess (gq—>"/q) dominates 7 production in all kinematic 
regions of pp scattering, as well as for low to moderate 
values of parton momentum fraction x in pp scattering. 
As mentioned previously, the gluon distribution is rela- 
tively well constrained at low x (x < 0.1) by DIS and 
DY data, but less so at larger x. Consequently, direct 
photon data have been incorporated in several modern 
global parton distribution function analyses and can, in 
principle, provide a major constraint on the gluon distri- 
bution at moderate to high x.c3'c3 

A pattern of deviations of direc±_photon data from 
NLO predictions has been observed, Ell however, with the 
deviations being particularly striking for the E706 exper- 
iment. c3 The origin of the deviations lies in the effects of 
initial state soft gluon radiation, or fc^. Direct evidence 
of this kr has long been evident from_T±Lfill-Yan, dipho- 
ton and heavy quark measurements. Ej'EHI The values of 
{kx) per parton vary from 1 GeV/c at fixed target ener- 
gies to 3-4 GeV/c at the Tevatron collider. The growth 
is approximately logarithmic with center of mass energy. 

Perturbative QCD corrections at the NLO level are 
insufficient to explain the size of the observed fcy and full 
resummation calculations are required to explain Drell- 
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Figure 14: The ratio of Rio with a jet cone size of 0.7 to Rio with 
a jet cone size of 0.4 for both CDF data and NLO QCD theory. 
Theory predictions are shown for several choices of scale. 



Figure 12: The ratio of (R 10 (Data) - R 10 {QCD)) to Rio (QCD) 
for NLO QCD predictions calculated using DYRAD. The parton 
distribution function set used is MRSA', and the baseline renor- 
malization and factorization scale is Mw Also shown are curves 
using other Q 2 scales. 
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Figure 13: The ratio Bio at Exmin = 30 GeV for both 0.4 and 
0.7 jet cones, compared to DYRAD predictions as a function of 
cts(Mz)- The CDF data are shown as horizontal bands. The the- 
oretical values for the MRSA' and CTEQ4 pdf families are shown 
as solid circles and open squares, respectively. 



Yan, W/Z and diphoton distributions. These resumma- 
tion calculations correctly describe the growth of the (kx) 
with center-of-mass energy. In a kx type of resummation 
calculation, there are typically two scales involved, which 
are important to the problem, and which differ greatly 
in magnitude from each other. In the case of direct pho- 
ton production, the two scales are the mass and px (or 
kx ) of the photon-jet system. Instead of examining the 
effects of the soft gluon resummation on the photon-jet 
mass cross section per se, one can instead look at the 
effects on the px distribution of the photon alone. 

Currently, there is no rigorous kx resummation cal- 
culation available for single photon production. The cal- 
culation is quite challenging in that the final state par- 
ton takes part in soft gluon emission and in color ex- 
change with the initial state partons, in contrast with 
the Drell-Yan and diphoton cases. Also, the calcula- 
tion is complicated by the fact that several overlapping 
power-suppressed contributions can contribute. In lieu 
of a rigorous calculation of the resummed direct photon 
px distribution, the effects of soft gluon radiation can be 
approximated by a convolution of the NLOj^xoss section 
with a Gaussian kx smearing function. C2'E3 The value 
of (kx) to be used for each kinematic regime is taken 
directly from relevant experimental observables, rather 
than from a theoretical prediction. 

The behavior of the kx smearing correction is quite 
different for the Tevatron collider and for fixed target 
energies. In Figure 15 is shown the comparison of NLO 



theory calculations (with and without the kx corrections) 
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Figure 15: The CDF and DO isolated direct photon cross sections, 
compared to NLO QCD theory without kx (dashed) and with kx 
enhancement for (kx) = 3.5 GeV/c (solid), as a function of px- 



Figure 16: The variation of kx enhancements (ratio of cross sec- 
tions with and without the kx corrections) relevant to E706 direct 
photon data at 31.6 GeV, for different values of average kx • Also 
shown is the kx correction for E706 used in the recent MRST fit. 



to the direct photon data from CDF and DO. The value of 
(kx) used (3.5 GeV/c) was takexL. directly from diphoton 
measurements at the Tevatron.Ej 

There are two points to note: (1) the agreement with 
the data is improved if the kx correction is taken into 
account and (2) the kx smearing effects fall off roughly 
as l/p T . The latter behavior is the expectation for such 
a power-suppressed type of effect. 

The kx correction obtained for E706 at a center of 
mass energy of 31.6 GeV is shown in Figure [jjj. The value 
of (kx) of 1.2 GeV was obtained from measurements of 
several kinematic variables in the experiment. cj The kx 
smearing effect is much larger here than observed at the 
collider and does not have the l/p T falloff. This can be 
understood from the following argument. At low px an 
(kx) value of 1.2 GeV/c is non-neglible compared to the 
px in the hard scattering, and the addition of the kx 
smearing both increases the size of the cross section and 
steepens the slope. At highpy (corresponding to large x), 
the unmodified NLO cross section becomes increasingly 
steep (due to the rapid falloff of the parton densities) and 
hence the effect of the smearing again becomes larger. 

For both the fixed target and collider cases, the kx 
smearing increases the size of the observed cross sec- 
tion. The direct photons are only measured above a given 
threshold. The direct photon px distribution is steep and 
the net effect of the kx smearing (or a more complete re- 



summation treatment) is to transport events from below 
the threshold to above the threshold, thus increasing the 
observed cross section. 

The uncertainty- in the value of (kx) is estimated by 
the E706 authors \B to be ±0.2 GeV. The effect of this 
variation on the kx correction is shown in Figure ^ and 
can be observed to be quite sizeable. Also shown is the kx 
correction used in the recent MRST paper. lJ The MRST 
paper uses an analytic kx smearing correction technique 
with an (kx) per parton value of 1.3 GeV/c. ^| 

The E706 direct photon cross sections for pBe col- 
lisions at y/s = 31.6 GeV is shown in Figure [l?] along 
with the NLJ^j-thcoretical predictions for the range of kx 
corrections. E§E3Ea Very good agreement is obtained with 
the use of the nominal value of (kx) ; the experimental 
cross section differs both in magnitude and shape from 
the NLO prediction with no kx correction. Also shown 
in the figure are the data and theoretical predictions for 
7r° production. NLO calculations for ir° production have 
a greater uncertainty than those for direct photon pro- 
duction since they involve parton fragmentation. How- 
ever, the kx effects are expected to be generally similar 
to those observed in direct photon production, and the 
7T° data can be used to extend tests of the consequences 



c The MRST paper quotes a smaller value for the (kx) used for 
E706. This smaller value is equivalent to 1.3 GeV/c in the conven- 
tion used in this talk and in Reference 29. 
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Figure 17: The photon and ir° cross sections from E706 compared 
to fc-p-enhanced NLO QCD calculations. In the bottom plot, the 
quantity (Data-Theory) /Theory is plotted, using fc^-enhanced cal- 
culations for several values of < kx >. The error bars have exper- 
imental statistical and systematic errors added in quadrature. 

of kx smearing. The presence of an additional kx value 
of a magnitude similar to that needed for single photon 
production leads to a substantially improved agreement 
between the ir° data and theory. 

The same comparison is made in Figure [l^ using 
the CTEQ4HJ pdf's. Good agreement is observed at 
low px but the theoretical prediction is larger than the 
data at high px- As mentioned previously, the dominant 
mechanism for direct photon production is gluon-quark 
scattering, so an increase in the gluon distribution in 
this range by a factor of 2 leads to an increase in the 
direct photon cross section by the same factor. Similar 
conclusions-.axfi obtained with the E706 data at y/s — 
38.8 GeV.M 

A decrease in the (kx) per parton at high px would 
lead to a better agreement of the E706 data with the 
CTEQ4HJ pdf predictions. There are several possible 
suppression mechanisms for soft gluon emission in this 
kinematic region that are not taken_into account in the 
simple kr model discussed above. cB (An experimental 
determination of the dependence of the < kx > as a func- 
tion of px/x is difficult due to the diminishing statistics 
at higher px-) This possibility is currently under inves- 
tigation. 

A comparison of the CTEQ4M and CTEQ4HJ gluon 
distributions and the gluon obtained by fitting the kx- 



Figure 18: The photon cross section from E706 compared to kx- 
enhanced NLO QCD calculations using the CTEQ4HJ parton dis- 
tribution set. 

corrected E706-data (along with the CTEQ4 DIS and 
DY data sets) c3 is shown in Figure |l9[ As might have 
been expected from Figure [l^, the gluon distribution ob- 
tained from this fit agrees well with the CTEQ4M gluon 
distribution and lies below the CTEQ4HJ gluon distri- 
bution at high x. Also shown in the figure is the gluon 
distribution from the MRST pdf, which incorporates the 
fc^-corrected data from WA70 in the fit. MThe larger 
kx correction for E706 observed in Figure [Lq implies a 
smaller gluon distribution as is observed in Figure |l9|. 

As has been discussed, there is a great deal of theo- 
retical uncertainty in the calculation of the fixed target 
direct photon cross sections. Even if the Gaussian kx 
smearing ansatz given above were formally correct, the 
uncertainty in the value of (kx) to be used leads to a 
large variation in the predicted cross section. This vari- 
ation makes the use of fixed target direct photon data in 
pdf fits somewhat problematic. A more complete resum- 
mation calculation should, with the appropriate experi- 
mental input, be able to predict the (kx) per parton for 
each kinematic condition and thus may be able to rescue 
the situation. 

There has also been much recent interest in study- 
ing the effects of resumming large logarithms of the form 

d ln the MRST paper, the E706 data are not used directly in the 
fits, but the fit results are compared to the-idata from E706 af- 
ter applying an appropriate kx correction. £3 Good agreement is 
obtained. 
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Figure 20: Diagrams indicating initial parton emission in ep scat- 
tering. 



be observed. As HERA has continued its successful op- 
eration, the data available for analysis by both of its ex- 
periments, HI and ZEUS, has increased steadily. Along 
with the increase in statistics has come an increase in the 
level of understanding of the detector systematics, allow- 
ing for more precise comparisons of data to perturbative 
QCD. 



Figure 19: A comparison of the CTEQ4M, MRST and CTEQ4HJ 
gluons and the gluon distribution derived from fits that use E706 
data. The gi and gl gluon densities correspond to the maximum 
variation in (kj<) that MRST allowed in their fits. 



ln(l — xt)-0 As xt approaches 1, for any hard scatter- 
ing process, the perturbative cross section is enhanced by 
powers of ln(l — xt) that have to be resummed at all or- 
ders. These types of effects should currently be negligible 
for direct photon production at the Tevatron collider (be- 
cause data are only available for relatively small values of 
xt), but may be important at fixed target energies. The 
net effect for E706Js a significant increase in the cross 
section at high pt-Ej A treatment that includes both kx 
and threshold resummation effects may be necessary for 
a more satisfactory description of the fixed target data. 
Recent, theoretical progress has been made in this direc- 
tion. c3 

In the CTEQ4 pdf fits (and the upcoming CTEQ5 
fits as well), the inclusive jet cross sections from CDF 
and DO provide an additional constraint on the gluon 
distribution at moderate to large x values. Because of the 
theoretical uncertainties mentioned above, the CTEQ5 
fits will not use fixed target direct photon data. 

4 HERA 

HERA is a positron-proton collider (27.5 GeV e + on 820 
GeV protons) with a total center of mass energy of about 
300 GeV. The large center-of-mass energy available at 
HERA offers a large phase space for the hadronic final 
state in DIS events, thus allowing clean jet structures to 



4-1 Parton Evolution Dynamics 

Of particular interest at HERA are measurements 
which discriminate aipong parton evolution schemes. 
DGLAPeS and BFKLcil describe the evolution towards 
large values of Q 2 and 1/xsj, respectively. DGLAP evo- 
lution resums terms of the form \n(Q 2 /Q 2 l ) and involves 
a strong ordering in kx of the gluon emissions with the 
hardest emissions occurring near the top of the gluon 
ladder. (See Figure pp|a.) In the BFKL model, terms 
of the form ln(l/a:) are resummed and gluon emissions 
are not ordered in transverse momentum fcT-i-4 solution 
of the parton evolution equation by CCFM e3 approxi- 
mates the BFKL equation in the low xb-j limit and the 
DGLAP equation in the high XBj limit. A lack of order- 
ing similar to that found in BFKL evolution has been in- 
corporated in the color dipole model. In the color dipole 
model, gluon emission originates from a color dipole that 
is stretched between the scattered quark and proton rem- 
nant. The result is a cascade of independently radiating 
dipoles with the gluons not ordered in kx- (See Fig- 
ure |2(i|c.) In addition, for processes where the photon 
structure is resolved, the hardest emissions given in the 
QCD matrix element may occur anywhere in the ladder, 
with increasingly soft emissions along the ladder towards 
both the proton and the photon. (See Figure 20 d.) The 
results presented previously for the Tevatron collider and 
fixed target measurements are all governed by DGLAP 
kinematics. At HERA, results have been obtained in all 
of the kinematic regimes discussed above. 

A wide variety of leading order (in a s ) theoretical 
programs are available for comparison to the HERA data 
incorporating the above evolution schemes: 
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Figure 21: The differential dijet cross sections from ZEUS at the 
parton level versus (a) Q 2 , (b) xgj, (c) £ and (d) V% ma ^l Q 2 ■ The 
inner(outer) error bars indicate the statistical error (statistical and 
systematic errors added in quadrature). The shaded area repre- 
sents the error due to the uncertainty on the jet energy scale of 
±3.5%. The solid line indicates the NLO QCD predictions from 
MEPJET. 



• BFKL calculations at t 
fragmentation functions 



parton level and with 



• the DGLAP-ba&ed parton shower Monte Carlo 
model LEPTO0 

• the CCFM-based Linked Dipole Chain model 
LDCB 

• the color dipole Monte Carlo model ARIADNE 
(which assumes a chain of independently radiating 
dipoles spanned by color-connected partonsLJ 

• the resolved photon model RAPGApEzI 

In addition, there are a number of DGLAP-based cal- 
culations in next-tctleading order of p., including DIS.- 
ENTE3, MEPJETEJ, DISASTER++E3 and JETVIPB, 
that are available in the form of flexible Monte Carlo 
generators. 

Comparisons between the data and the theory can 
be made either at the detector level (with an appropri- 
ate Monte Carlo simulation of the detector response), the 
hadron level (with corrections for acceptance and migra- 
tion) or at the parton level (correcting for hadronization 
effects). 
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Figure 22: The dijet cross sections from HI double differential in 
Q 2 and £. Also shown are the NLO QCD predictions of DISENT. 



4-2 Dijet Cross Sections in DIS 

In DIS (in the naive quark-parton model), the virtual 
photon is absorbed by a single quark or antiquark in 
the proton. This results in one jet from the struck 
quark(antiquark) and one jet from the proton remnant 
(1+1 configuration). To first order in a s (lowest order 
for dijet production), two jets (in addition to the proton 
remnant jet) with balanced transverse momentum are 
produced in the photon-proton center-of-mass(2+l con- 
figuration). There are two subprocesses responsible for 
dijet production in DIS: boson-gluon fusion (j*g—>qq) 
and QCD Compton scattering (7*9— ><?<?)■ At low x and 
Q 2 , the large gluon density leads to the dominance of 
the boson-gluon fusion subprocess and allows for a di- 
rect sensitivity for the gluon distribution in an x region 
below the fixed target direct photon experiments. The 
presence of the strong interaction vertex may allow a 
measurement of a s . Jet measurements have been con- 
ducted with variations of both the iterative cone and fcy 
algorithms, and in the lab, Breit and hadron center-of- 
mass reference frames. 

ZEUS has submitted to this conference a measure- 
ment of dijet production in DIS events from a data 
sampk^CQEiesponding to an integrated luminosity of 6.4 
pfc-^BEaQ Cuts were placed on Q 2 (7-100 GeV 2 ), y(> 
0.04) and the scattered positron energy (> 10 GeV). The 
jets were reconstructed with a cone algorithm (R=l) and 
were required to have a transverse energy greater than 4 
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Figure 23: The dijet cross sections from HI double differential in Q 2 
and xgj . Also shown are the NLO QCD predictions of D1SENT. 



GeV/c in both the laboratory and center-of-mass frames, 
and to have a pseudorapidity in the laboratory frame in 
the range from -2 to +2. The differential cross sections, 
corrected to the parton level, are shown in Figure 
The corrections applied are typically 20-40%. The vari- 
able £(= x_bj(1 + m'jj/Q 2 ) is of particular interest since 
it is related to the momentum fraction of the quark or 
gluon emitted from the proton in the leading order QCD 
picture. 



As can be observed in Figure 21, the data span the 
range from 0.2 < p\jQ 2 < 30, so the events are char- 
acterized by two scales, p\ (square of jet transverse mo- 
mentum) and Q 2 . A comparison of an exact NLO QCD 
calculation (MEPJET) to the data demonstrates that the 
theory adequately describes the shape of all of the dis- 
tributions but has a normalization off by approximately 
34%. The difference is comparable in size and of the same 
sign as the hadronization correction applied to convert 
the data to the parton level. 

The difference between theory and experiment ob- 
served by ZEUS might originate in the effects of soft 
gluon radiation in the regime (same minimum transverse 
momentum cut on both of the jets) where such effects 
may be important. It has been proposed to measure the 
dijet rate using asymmetric cuts on the transverse mo- 
menta of the two jets in order to minimize this type of 
correction to the NLO calculation. E3 Placing an addi- 
tional cut on the larger of the transverse momenta of the 



two jets decreases the magnitude of both the predicted 
and measured dijet cross sections, but the relative differ- 
ence between the two still remains substantial. Similarly, 
increasing the cuts on the transverse momenta of the two 
jets results in both the measured and calculated cross 
sections decreasing, withJjhe ratio of the two remaining 
approximately constant, til 

HI has measured the dijet cross section in DIS events 
at HERA and from this measurement has extracted a 
determination oilhc^gluon distribution in the x range 
from .01 to 0.1.E3H'Ej The analysis involves a large data 
sample (36 pb" 1 ) and utilizes fcx jet algorithms for jet de- 
termination in the Breit frame. The Breit frame, where 
the virtual photon collides head-on with the scattered 
quark, is well-suited for studies of dijet production. In 
this frame, the jet transverse energy directly reflects the 
hardness of the underlying QCD process. Three differ- 
ent variations of the kr algorithm were utilized with 
the smallest hadronization corrections being present in 
the longitudinally invariant Ut algorithm. The double- 
differential dijet cross sections (corrected to the hadron 
level) are shown in Figure ^2] as a function of Q 2 and 
£ and in Figure ^3| as a function of Q 2 and XBj, com- 
pared to the NLO program DISENT. Good agreement 
is observed but the data have not been corrected for 
hadronization effects (which, however, are expected to 
be small). 

Most of the data presented above by HI are at higher 
Q 2 than the previously discussed ZEUS measurement. 
There is some overlap, though, and the reason for the 
agreement by HI with NLO theory and the disagreement 
by ZEUS is still under study. It may be related to the two 
studies being carried out in different frames of reference, 
with different jet algorithms and the jet corrections being 
to the hadron level for HI compared to the parton level 
for ZEUS. 

To determine the gluon distribution, the dijet cross 
section distributions d 2 a / dQ 2 d£ and d 2 a/dQ 2 dxBj were 
utilized (in the range 200 < Q 2 < 5000 GeV 2 , where 
the NLO theory describes the data well) along with the 
HI neutral current DIS data (to determine the quark 
densities in this x range). In the fit, the data are com- 
pared to the product of the NLO QCD prediction and 
the hadronization correction. The results of the gluon 
fit are shown in Figure |24| for the x range 0.01 to 0.1 
with a factorization scale (/x 2 ) of 200 GeV 2 . The result 
is slightly higher than the gluon distributions obtained 
from several global analyses (although compatible within 
the error) and is in good agreement at x = 0.01 with the 
gluon obtained from the QCD analysis of the HI F2 data. 

HI has also reported a determination Q_Lthe strong 
coupling constant a s from the dijet sample. E3E3 A mod- 
ified JADE kr jet algorithm was applied to a sample of 
NC DIS events in the Q 2 range from 200- 10000 GeV 2 . 
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Figure 24: The error band of the gluon density in the proton from a 
NLO QCD fit to the HI dijet cross sections. The result is compared 
to the gluon densities from different parton distribution sets and 
the result from a fit to the HI structure function data. 



Figure 25: The distribution of the differential jet rate 1/2 corrected 
for detector and hadronization effects compared to the NLO predic- 
tion from MEPJET. The full line shows the NLO prediction using 
the fitted value of a s . 



The restriction of Q 2 > 200 GeV 2 provides better ac- 
ceptance for the final state jets and restricts the range of 
initial parton x to large values where the parton densities 
are better known. 

The jet algorithm calculates the scaled quantities 
mfj/W 2 of pairs of calorimeter clusters (i,j), where W 2 
is the total invariant mass of all clusters and my is the 
mass of clusters i and j. The clusters with minimum 
rriij/W 2 are added together; this procedure is repeated 
until exactly (2+1) jets remain. The smallest scaled jet 
mass given by any combination of the (2+1) jets is de- 
fined to be the observable 3/2- A cut on y-i (yi > 0.01) is 
imposed to increase the fraction of events with a clear 
(2+1) jet structure, thus enhancing the sensitivity to 
a s . The HI yi distribution, normalized to the num- 
ber of DIS events in the kinematic sample and corrected 
for detector and hadronization effects, is shown in Fig- 
ure H[ along with the predictions of the NLO program 
MEPJET. Predictions are shown for A^^L = 100 MeV 

A/ S 

and 600 MeV, along with the best fit value of 320 MeV. 
This value of A corresponds to a value of a s (M§) of 
0.118 ± Om2to°%(systematic). There is an additional 
theory systematic error of to'oog- This result is in agree- 
ment with the world average, albeit with a large error. 
The largest experimental systematic error is due to the 
model dependence of the detector and hadronization cor- 
rections; one of the largest sources of theory uncertainty 
is due to the imperfect knowledge of the pdf's (and in 
particular the gluon distribution) in this kinematic range. 
(Gluon-initiated processes account for about 50% of the 



(2+1) jet events used in the analysis.) There is a strong 
correlation between the fit value of a s and the size of the 
gluon distribution as was observed in determinations of 
a s in jet production at the Tevatron.liJH 

4-. 3 Forward Jet Production at HERA 

One of the significant discoveries made at HERA was the 
steep rise of the proton structure function F2(x,Q 2 ) in 
the region of small x (x < 10 -3 ). In the BFKL approach, 
the leading terms in ln(l/a:) which appear together with 
the In Q 2 terms in the evolution equation are resummed. 
The BFKL terms may lead to a steeper F2 behavior, 
but from the existing F2 data, it is not possible to un- 
ambiguously determine whether the BFKL mechanism 
plays a role in the HERA x range. The BKFL mech- 
anism predicts additional contributions to the hadronic 
final state from high transverse momentum partons trav- 
elling forward in the HERA frame. These forward-going 
partons may be detected experimentally as jets and may 
result in an enhancement of the forward jet cross section 
when compared to cither exact NLO QCD calculations 
or parton shower model calculations based on DGLAP 
evolution. This is the same type of BFKL physics that 
was discussed earlier for large dijet rapidity separation at 
the Tevatron. In this case, the large rapidity separation 
is between the current and forward jets. 

As previously, there are two hard scales relevant for 
low x forward jet production: the squared momentum 
transfer of the photon, Q 2 , and the squared transverse 
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Figure 26: The hadron level forward jet cross section as a function 
of E^,j et /Q 2 from ZEUS. The data are compared to the RAP- 
GAP Monte Carlo model with direct and resolved contributions, 
to LEPTO and to ARIADNE. The shaded band corresponds to the 
uncertainty from the energy scale of the calorimeter. 



energy of the jet E\- et . The variable E^- et /Q 2 can be 
varied from very small values to very large values. A 
low value of E\jQ 2 corresponds to the standard DIS 
regime where DGLAP dynamics is dominant. In the 
range where E^/Q 2 is approximately 1, BFKL dynam- 
ics becomes important (and DGLAP parton evolution 
is suppressed) ^] while in the regime where E\ > Q 2 , 
the jet begins to probe the structure of the photon. 
The forward jet cross section from ZEUS— (corrected to 
the hadron level) is shown in Figure |2(| E3'Ea All three 
Monte Carlo programs shown describe the data well for 



E 2 /Q 2 « 1. ARIADNE and RAPGAP work in the 
BFKL region {E\ = Q 2 ) while only RAPGAP is suc- 
cessful in the regime where E\ » Q 2 . The RAPGAP 
Monte Carlo model contains resolved as well as direct 
photon contributions. A resolved virtual photon con- 
tribution could account for the excess of forward jets 
with respect to the standard DGLAP models. In Fig- 
ure the hadron level forward cross section is plotted 
as a function of Bjorken x. The agreement with RAP- 
GAP (with both resolved and direct components) is ex- 
cellent. However, the amount of resolved contribution 
to the forward jet cross section has a wide range of un- 
certainty that makes definitive conclusions difficult. The 
shaded band in Figure |2?] indicates the variation in the 
RAPGAP prediction when the factorizaton scale is var- 
ied from ^t 2 = E 2 jet /2 + Q 2 to /i 2 = AE 2 jet + Q 2 . Due to 
the large scale dependence, a comparison to exact NLO 
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e Specifically, BFKL dynamics is important in the forward jet region 
when E^/Q 2 is near one and Xj e t >> %Bj- 



Figure 27: The hadron level forward jet cross section from ZEUS as 
a function of xgj . The data are compared to the RAPGAP Monte 
Carlo model with direct and resolved contributions and to LEPTO. 
The shaded band on the top indicates the uncertainty due to fac- 
torization scale variation for the full RAPGAP prediction while the 
shaded band on the bottom indicates the same uncertainty for the 
direct contribution alone. 



calculations is needed but the appropriate hadronization 
corrections are large. 

HI has measured forward jet and dijet production 
requiring a transverse-energy larger than 3.5 GeV/c in a 
cone of radius 1.0. BEJ A cut of 0.5 < E 2 jet /Q 2 < 2 is 
applied to enhance BFKL effects. The forward jet cross 
section is shown in Figure ^8], as a function of x compared 
to several models. ARIADNE and RAPGAP (with a re- 
solved photon contribution) lie closer to the data than 
does the DGLAP-based program LEPTO. NLO parton 
calculations using DGLAP parton densities, as for exam- 
ple DISENT, disagree with the data both in shape and 
normalization. The numerical BFKL calculations at the 
parton level lie above the data but describe the shape 
fairly well. 

The forward dijet cross section is measured by HI to 
be 6.0±0.8(stat)±3.2(sys) pb, in agreement with the pre- 
dictions from ARIADNE and RAPGAP. The dijet cross 
section is roughly 1% of_the forward jet cross section. Re- 
cent BFKL predictionsEil predict on the other hand that 
3-6% of the forward jet data should contain 2 or more 
forward jets. 

Single high px particle production can also be used as 
a probe of QCD dynamics in the forward region. The am- 
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Figure 28: A comparison of the HI forward jet cross section to 
various QCD calculations. 
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Figure 30: The measured three-jet cross section with respect to the 
three-jet invariant mass. The inner error bar gives the statistical 
error and the outer the sum of statistical and systematic errors in 
quadrature. 




Figure 29: The forward tt° cross sections as a function of Xgj for 
four different values of the variable = / E pro ton Here n n is the 
number of 7r°'s and N is the number of events in the distribution. 
Comparisons of four different QCD models are overlayed. For the 
plot on the lower right, a BFKL calculation is shown and has been 
been divided by two for presentation. 



biguity of jet definition implicit in the forward jet search 
is missing, and the smaller spatial extent allows the mea- 
surement to smaller angles. In Figure |2^, the forward tt° 
data from HI is shown plotted vs Xgj for several values 
of Xt T {=E 1t I 'Ep roton ). The observed rise with decreasing 
XBj again provides evidence of more hard partonic radia- 
tion than predicted by DGLAP type Monte Carlo models 
(such as LEPTO), and is more reasonably represented by 
RAPGAP, where the photon acts as a resolved object. 

4-4 Multijet Photoproduction at HERA 

The study of multijet photoproduction provides a direct 
test of perturbative QCD predictions beyond leading or- 
der. Multijet kinematic observables have been previously 
studied for 3-6 jet production at FermilabE3 and for 2 jet 
production at HERA. El ZEUS has now measured, .3. jet 
final states in photoproduction events at HERA.Ej'Ej 

Photoproduction events were selected by restricting 
the transverse momentum to the positron to be less than 
1 GeV and the photon-proton center-of-mass energy to 
be in the range from 134 to 269 GeV. Jets were defined 
using a kx cluster jet algorithm with the first two jets 
having a transverse energy greater than 6 GeV/c and 
the third a transverse energy greater than 5 GeV/c. The 
requirement of relatively high transverse energy for the 
jets ensures that the process can be calculated by per- 
turbative QCD. 

The three-jet invariant mass distribution is shown 
in Figure BO and compared to the order aa 2 s calcula- 
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Figure 31: The distributions of the angles 005(63) and 1/13. The 
thick lines show the fixed order perturbative QCD calculations and 
the thin lines represent the parton shower Monte Carlo predictions. 
The dotted curve shows the distribution for a constant matrix ele- 
ment. 
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Figure 33: The inclusive charged particle distribution from ZEUS, 
in the current fragmentation region of the Breit frame. The inner 
error bar is the statistical and the outer error bar shows the statis- 
tical and systematic errors added in quadrature. The open points 
represent data from e+e - experiments divided by two to take into 
account q and q production. 



Figure 32: The measured cross section do /dips and the normalized 
tp 3 distributions for ZEUS. The PYTHIA and HERWIG predictions 
are shown by solid and dashed lines respectively and the prediction 
from PYTHIA with color coherence switched off is shown by the 
dashed-dotted line. 

tions. The two calculations are leading order for the 
variable under study (since there are three jets in the fi- 
nal state) but still provide good agreement with the data. 
The Monte Carlo programs PYTHIA and HERWIG con- 
tain only 2 — > 2 matrix elements but a third jet can be 
provided by gluon radiation. The Monte Carlo programs 
predict the correct shape for the cross section but have 
a normalization too low by 20-40%. 

The distributions for 63 (the angle between the high- 
est energy jet and the beam direction) and 1/73 (the angle 
between the plane containing the three jets and the plane 
containing the highest energy jet and the beam direction) 
are shown in Figure 

The cos (#3 ) distribution shows the expected Ruther- 
ford scattering form ((1 — cos^)) -2 ). Both the #3 and 
1/73 distributions differ dramatically from phase space and 
agree well with both the fixed order calculations and the 
Monte Carlo models. 

The QCD phenomena of color coherence can be 
tested using the '03 distribution. In Figure ^2| is shown 
the ip3 data again, along with predictions of PYTHIA 
and HERWIG (both with color coherence implemented) 
and PYTHIA (with color coherence turned off). Color 



coherence disfavors gluon radiation into certain angular 
regions which are determined by the color flow of the pri- 
mary scatter. PYTHIA with no color coherence predicts 
a much flatter ^3 distribution than observed in either the 
data or in HERWIG and (the color coherence version of) 
PYTHIA. 

4-5 Jet Fragmentation at HERA 

Fragmentation functions characterize the process of 
hadron formation in jet production and decay. A nat- 
ural frame to examine the details of jet fragmentation in 
DIS events is the Breit frame, defined previously. The 
current region in the Breit frame is analagous to a sin- 
gle hemisphere in e + e~ collisions and the fragmentation 
properties of these quarks can be directly compared to 
the fragmentation of the struck quark in the proton. The 
ep Breit frame equivalent of the e + e~ scaled hadron mo- 
mentum, x p = 2p hadron / y /s 1 is x p = 2p hadron /Q, where 
only hadrons in the current hemisphere are considered. 
The fragmentation function from the ZEUS experiment, 
plotted as a functiorLolQ for different intervals of 
shown in Figure |33|. P 8 ^ 9 ! Approximate scaling is observed 
at moderate values of x p with clear scaling violations be- 
ing apparent at both large x p (decrease with Q 2 ) and 
small Xp (increase with Q 2 ). The HERA data overlap 
the kinematic range of the e + e~ data and good agree- 
ment between both types of experiments is observed. 
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Figure 34: The charged particle distributions l/Ndn/dln(l/x p ) as 
a function of Q. Only statistical errors are shown. 

The small x p region is better investigated using the 
variable £ = ln(l/iEp). The Modified Leading Logarithm 
Approximation (MLLA) together with Local Parton- 
Hadron Duality (LPHD) predict both the shape of the £ 
distribution ("hump — backed") and the evolution of the 
peak and the width of the £ distribution with energy. 119 
The charged particle distributions are plotted as a func- 
tion of £ in Figure |34| for various Q values, and the values 
of Cpeafe and £, W idth are shown as a functkin of Q(a/s) for 
HI (e + e~ experiments) in Figure |35|. EaO The observed 
peak and width of the £ distributions at HERA agree well 
with the e + e _ data and with the MLLA predictions. 

4-6 Jet Shapes 

For cone jet algorithms, a useful representation of the 
internal structure of a jet is given by the jet shape. At 
sufficiently high energies, the jet shape should be calcu- 
lable in perturbative QCD, with gluon jets broader than 
quark jets. At HERA, jet production has been observed 
in both neutralpiiuricnt (NC) and charged current (CC) 
DIS at high Q 2 .Ea€2l As mentioned previously, in DIS, jet 
studies can be carried out in different frames. The appro- 
priate frame for jet shape studies is still under discussion 
and the effects of boosting to different frames have not 
been fully investigated yet. 

ZEUS has carried out a comparison of jet shapes 
in NC and CC interactions at Q 2 > 100 GeU 2 _along 
with jets from e + e~ and pp and "fp collisions. O Jets 



Figure 35: A comparison of HI results showing the evolution of 
the peak (a) and the width (b) of the fragmentation function as a 
function of Q. Also shown are e + e — results at the corresponding 
values of center-of-mass energy. The solid line is a fit of HI data 
alone to MLLA/LPHD expectations. 

are measured with an iterative cone algorithm, in the 
laboratory frame of reference, with a radius R of 1 and 
an Et value of greater than 14 GeV. The jet shape was 
measured with the ZEUS calorimeter and corrected to 
the hadron level. 

The integrated jet shape ip(r) (the fraction of jet 
energy inside a cone of radius r compared to the total 
inside radius R) is shown in Figure |36| for NC events, re- 
solved photoproduction events (x° bs <_p.75) and direct 
photoproduction events (x° bs > 0.75).O In direct photo- 
production events, the photon acts as a point-like object 
while in resolved events the parton structure of the pho- 
ton is probed. I* measures the fraction of the photon's 
energy that goes into the two highest Et jets. The jets 
produced in NC DIS are narrower than those in dijet 
photoproduction, but closer to those dominated by di- 
rect processes. In NC DIS events, most of the final-state 
jets are quark jets (e + q— >e + q); direct photoproduction is 
dominated by the subprocess jg^qq but has contribu- 
tions from the subprocess ~/q—>qg- The resolved photo- 
production events have a larger fraction still of final-state 
gluon jets as evidenced by the larger jet width. 

Jets from NC and CC DIS with Q 2 > lOOGeU 2 and 
jet Et values in the range from 37-45 jQeV ace compared 
to jets of similar Et values from CDFO, D0L3 and (from 
e + e~ collisions) OPALL3 in Figure |3?]. For all three ex- 
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Figure 36: The measured integrated jet shape corrected to the 
hadron level from ZEUS in NC DIS, and in resolved and direct 
photoproduction. 



Figure 37: The measured integrated jet shapes corrected to the 
hadron level in NC and CC DIS events at ZEUS and in pp collisions 
at CDF and DO and from e + e - collisions at OPAL. 



periments, an iterative cone algorithm with a radius R 
of 1 is used. For the two collider experiments, an under- 
lying event level corresponding to minimum bias events 
is subtracted. No such correction is needed for the e + e~ 
events, or for the HERA DIS events (in the kinematic 
region being considered). The jets from the Tevatron 
Collider are significantly broader than the DIS jets from 
HERA and the e + e~ jets from OPAL. This difference 
is primarily due to the larger fraction of gluon jets ex- 
pected at the Tevatron with perhaps some small part of 
the difference being due to the extra contributions to the 
underlying event that may be present at the Tevatron. 

HI has carried out a measurement of the internal jet 
structure in an inclusive DIS dijet sample in the kine- 
matic domain UL < Q 2 < 120 GeV 2 and 2xl0~ 4 < 
x Bj < 8xlCT 3 . El Jets were reconstructed in the Breit 
frame using both the iterative cone and kr jet algo- 
rithms, with the requirement that ETBreit > 5 GeV/c 
and —1 < rjj et ,iab < 2. Jets become more collimatcd 
as ErBreit increases, with the dependence becoming 
more pronounced for the cone algorithm. For constant 
ErBreit, jets are narrower (broader) towards the photon 
(proton) hemisphere. These dependences become smaller 
as ETBreit increases. A possible explanation for this be- 
havior is that the internal jet structure is influenced by 
particles close to or produced by QCD radiation near the 
proton remnant. 

Jets defined by the kx algorithm tend to be more 
collimatcd than those defined by the cone algorithm. 
The dependence on ETBreit and ijBreit is also stronger 



for the cone algorithm. Generally, the jet shapes are 
well described by QCD Monte Carlo predictions with 
LEPTO having a tendency to produce broader jets to- 
wards the proton remnant direction, HERWIG producing 
jets which are too narrow (especially at large Et Breit and 
V Breit), and ARIADNE lying between the two and giving 
a good overall description of the data. 

5 Conclusions 

DGLAP-based perturbative QCD calculations have been 
very successful in describing data involving jets and pho- 
tons at both the Tevatron and at HERA. Most of the ar- 
eas in which the remaining disagreements/controversies 
exist involve either an uncertainty in the gluon distri- 
bution (the high Et jet cross section at the Tevatron) 
or the influence of two scales in the measurement (fixed 
target direct photon production and forward jet produc- 
tion at HERA). For the case of forward jet production, 
BFKL effects may or may not be important; a proper 
treatment of photon structure seems to describe the data 
both the BFKL region and beyond. Recent calculations 
of the forward jet production cross sections based on the 
BFKL approach showed unusually large next-to-leading 
order corrections, raising the question on their predictive 
power. L3 A deeper understanding of the origin of these 
large corrections is needed before a comparison to the 
data may be meaningful. 

For fixed target direct photon production (particu- 
larly in the case of E706), soft gluon effects are extremely 
important, changing both the shape and normalization 
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of the cross section with both the kx and Sudakov re- 
summation formalisms being required. There has been 
a great deal of theoretical effort on this problem; a suc- 
cessful resolution will allow the quantitative treatment of 
direct photon data in pdf fits again, provide a window on 
an area of very interesting physics, and finally, settle the 
question of the large x gluon distribution. 

Both CDF and DO are undergoing major upgrades 
for Run 2, which is scheduled to begin in the spring of 
2000. Each experiment will put in place a greatly im- 
proved detector and will accumulate (from 2000 to 2003) 
a data sample on the order of 2 a factor of 20 in- 

crease over Run 1. A data sample this size will enable 
the high Et jet cross section to be probed in much more 
detail as well as allowing a variety of Tevatron QCD mea- 
surements to be performed with greater precision. 

HI and ZEUS will continue the analysis of the data 
taken with positrons in 1994-1997. HERA switched to 
electron running this year and plans to deliver approx- 
imately 60 pb^ 1 over the next two years. In 2000, the 
HERA machine will be upgraded for high luminosity run- 
ning, with yearly rates of 150 pb^ 1 expected. 
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